Entanglement, combinatorics and finite-size effects in spin-chains 
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We carry out a systematic study of the exact block entanglement in XXZ spin-chain at A = —1/2. 
We present, the first analytic expressions for reduced density matrices of n spins in a chain of 
length L (for n < 6 and arbitrary but odd L) of a truly interacting model. The entanglement 
entropy, the moments of the reduced density matrix, and its spectrum are then easily derived. We 
explicitely construct the "entanglement Hamiltonian" as the logarithm of this matrix. Exploiting the 
degeneracy of the ground-state, we find the scaling behavior of entanglement of the zero-temperature 
mixed state. 



Introduction. Entanglement is a central concept in 
quantum information science and it is becoming a com- 
mon tool to study and analyze extended quantum sys- 
tems because of its use in detecting the scaling behav- 
ior close to a quantum critical point [ij. It is has been 
pointed out that this scaling behavior is connected with 
the efficiency of numerical methods as quantum and den- 
sity matrix renormalization group (DMRG) [l[ . 

Let p be the density matrix of a system and let the 
Hilbert space be written as a direct product H = Ha "X) 
Ti-B- ^'s reduced density matrix is pA — Tr^p and the 
entanglement entropy is the corresponding Von Neumann 
entropy 



Sa = -Tr APA^og PA , 



(1) 



and analogously for Sb- When p corresponds to a pure 
quantum state Sa — Sb- 

When A is a segment of length n of an infinite one- 
dimensional system in a critical ground state, the corre- 
sponding entanglement entropy iSyidiverges as the loga- 
rithm of the sub-system size 0, S [3| 



log n 



(2) 



where c is the central charge of the associated confor- 
mal field theory (CFT) and c'l a non-universal constant. 
Away from the critical point, Sn saturates to a constant 
3 1 proportional to the logarithm of the correlation length 
J|. These properties made the entanglement entropy a 
basic tool to analyze ID models. While it is impossible 
to mention here all the important contributions in the 
field, we refer the interested reader to the reviews 

In recent times, it has been remarked by few authors 
m that the reduced density matrix p„ contains much 
more information than Sn- To the best of our knowledge 
the full reduced density matrix is known only for free 
systems Q and is difficult to be obtained by numerical 
methods because these tend to focus on its eigenvalues. 
In order to go beyond free systems and to study the ef- 
fect of strong interactions, in this letter we report a first 
systematic study of the reduced density matrix of the an- 
tiferromagnetic XXZ chain at A = —1/2 defined by the 



Hamiltonian 



(3) 



with periodic boundary conditions (cTi+i = CTi) and an 
odd number of sites. Here a^'^'^ stand for the Pauli ma- 
trices at the site i. This critical model has the unique 
property that all the components of the ground-state 
wavefunction are integer multiples of the smallest one 
We will argue in the following that this property 
suffices to get pn{L) for n < 6 and arbitrary L: an ex- 
ceptional result for a truly interacting system. Another 
unique feature of this spin-chain is that the ground-state 
energy, which is doubly degenerate for any finite L has no 
corrections to scaling: ~ — (3/2)L exactly. We refer 
to the two ground states as \^±) with the upper index 
being the sign of the total spin in the z-direction. 

This work is complementary to two recent papers. In 
Ref. pn has been calculated in the thermodynamic 
(TD) limit (also for n up to 6). In Ref. [l3| the connec- 
tion with loop-models has been explored, for a different 
measure of entanglement. In contrast, here the emphasis 
is on combinatorial and finite-size scaling aspects. The 
exact value of some elements of the reduced density ma- 
trix for smaller values of n is also known for general A 
[m and all of them up to n = 6 for A = —1 [l2 |. 

Analytic results for p„. Since the ground-state en- 
ergy is exactly proportional to the system size, and the 
Hamiltonian is represented by a matrix with rational el- 
ements, also the ground-state vector has only rational 
components. Suitably normalized, all ground-state com- 
ponents are integer multiples of the smallest one. The 
ground-state for system sizes up to L = 25 can be then 
obtained with absolute precision in very modest com- 
puter time. With these ground-states we can construct 
the corresponding density matrices Pn{L) with n < L. 
Any element of Pn{L) is necessarily a rational number 
and in fact a rational function of L, with numerator 



72]. 



and denominator of degree 

guess the denominator to be 2"^L" nfc=V^ {L'^—Ak'^) 
(Anti)symmetrized with respect to the two ground states 



The data suffice to 

L"/2J/r2 Au2\n-2k 



2 



Pn{L) turns out to be an even (odd) function of L. As a 
result it can be determined completely for general L and 
for n < 6. For example, for n = I and Ji — 2 we obtained 



L + 1 
L-1 



p,{L) = ^ X 

/2(L + 2)2-2 

6L^-6 + 3 

5L2 + 3 6L2-6 



(4) 
(5) 
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2(L-2)2 -2 / 



The other reduced density matrices are too large to dis- 
play. We enclose an electronic Mathematica file (the den- 
sity matrix is rho[L,n] with L odd and n — 1, • • ■ ,6). 
In the mathematica file and in the following, the indeces 
of the matrix are the decimal form of the binary number 
representing the site product state (1 for -f and for — ). 

We observed the following properties of p„: 
(i) The first and the last elements correspond to the prob- 
ability of a string of equal spins, i.e. the emptiness for- 
mation probabilities (EFP), E±{L,n), which is 



n 

k=0 



k\ {2,k + l)\ {L-k~-l)\ (^-f fc)! 
(2fc)! (2fc + 1)! (L + k)\ (-^ - fc - 1) 



(6) 



for minority (see [8|]) and majority spins respectively. 
These two elements approach the same limit as L — > oo, 
and are equal to p[l, 1] = An/2" , where An is the num- 
ber of n X n alternating sign matrices (ASM). 

(ii) Some elements satisiy relations that connect p„ to 
Pri-i when summing over one spin (see the appendix 
A of Ref. (§1). These equations can be used to de- 
rive expressions for some more elements. For exam- 
ple, we have /9„+i[l, 1] + p„+i[2, 2] = p„[l, 1], that com- 
bined with the ASM sequence for 1], gives p„[2, 2] = 
(2"yl„_i - A„)/2" . Similar relations can be derived for 
a few other elements. 

(iii) For L — > cxi all non-zero elements remain non-zero 
and reproduce the results of Ref. [i| . 

(iv) The analytic continuation to general L satisfies 
Pn{L)[i,3] = p„(-L)[2" + 1 - 2" -f 1 - j]. 

The entanglement entropy in the TD limit is 

S"! = log2, S'2 = 0.950749, 5*3 = 1.09287, 

54 = 1.19076, = 1.26588, Se = 1.32701 (7) 

the same as in Ref. Q. We are in position to study 
the finite-size effects. In Fig. [T] we plot Sn{L) — 
l/31ogL/7r, and we compare it with the CFT prediction 
[J| 1/3 logsin(7rn/L) -|- c[ valid for large enough n. We 
notice that all the results fall on the same curve, except 
for small deviations at rt = 1 . It is impressive and maybe 
unexpected that the asympotic scaling sets in for such 
small value of n. 




0.4 0.6 
n/L 

FIG. 1: Finite size scaling of the entanglement entropy S„{L) 
against the CFT prediction 1/3 log sin(7rn/I/) + c'j (full lino). 
We fixed c'l = 0.7305 jlj. 



Combinatorics. It has been noted in the study of 
Hamiltonian ([3]) that many physical quantities are se- 
quences of integers or rationals that can be recognized 
in terms of known ones, with the ASM sequence for the 
ground-state components that is the best known but 
not the only one [8, 15, ijj . 

Once a sequence has been recognized, the critical ex- 
ponents can be derived exactly from its asymptotic be- 
havior, which in part motivates the great interest in this 
kind of studies. One could have hoped that the eigen- 
values of the reduced density matrices are rational or at 
least simple functions of L, but this is not the case. It 
is then natural to move our attention to the elements 
of the density matrix itself. Because of the conserva- 
tion of the spin in the z direction, the density matrix 
has non zero entries only between states with the same 
magnetization {— ~n/2, —n/2 + 1, . . . n/2 ~ 1, n/2) and 
organizes into sectors, as evident from the block struc- 
ture above. We already reported the formula for the first 
and last element of pn{L), that are the only elements in 
the sectors with spin ±n/2. For the other sectors, the 
elements grow too quickly with n to recognize any sim- 
ple behavior. For this reason we explored the possibility 
that some non-trivial combinations of them could have a 
simpler structure. Let us start by considering the sectors 
with total spin ±(n/2 — 1). After some tries, we found 
that reasonable growing sequences are given by the trace 
of the matrix multiplied with the matrix (—!)*+■'. In the 
TD limit, this results in the following sequence: 



1, 2, 11, 72, 806, 14352, 



(8) 



for both sectors, up to an overall factor 2"" . This se- 
quence grows with n mildly, but we have been unable to 
recognize it. At finite L, the same trace for Pn{2n + 1) 
gives 

1, 2, 31, 124, 489, 1826, 6843, 25712, 97213, 

369478, 1410831, 5408272, (9) 
and 2, 2, 2, 2, 2, 2, 2, 2, • • • (10) 
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for the two sectors with Sz = ±(n/2 — 1) of Pn{'2n + 1) 
respectively, up to the overaU factor E^{2n + l,n) from 
Eq. The first sequence also grows in a reasonable 

way, but still we have not been able to guess it. For sec- 
ond one, the recognition is instead obvious and provide 
a non-trivial relation whose physical origin is still un- 
known. Similar sequences are easily derived for the other 
sectors. We want to stress here that this generation of 
sequences is not only an academic game. If we would 
have been able to find a sufficient number combinations 
of elements of the reduced density matrix that can be 
recognized, we would have been access to analytic forms 
of the reduced density matrix for any n and L. By simply 
looking at the elements of the density matrix this could 
seem an impossible task, but we explicitly showed that 
at least one combination of them is very easy and that 
there are other sequences that do not look prohibitive. 
The main reason why we reported these series here is 
to stimulate further studies in this direction in order to 
eventually achieve the knowledge of the full density ma- 
trix. 

The moments of Pn{L) 



i?(r)(L) = Trp«(L), 



(11) 



for a integer are sequences of rationals. For critical sys- 
tems they display the universal asymptotic behavior 0, 



L 7rn 
— sm — 

TT L 



- — c(q — 1/q)/6 



(12) 



from which one can reconstruct the full spectrum of 
Pn (L) 6] and its behavior is connected with the accuracy 
of some numerical algorithms based on matrix product 
states ^ . Despite this universal behavior has been 
only marginally considered in the literature 0, . In the 
TD limit, since all elements of /9„(oo) have a common de- 

2 

nominator 2" , all moments can be written as 



ri") 2-""' , with ri") inte; 



gers. 



(13) 



For example, the values of ri^' up to n = 6 are 



(2) 
„(2) 



130, 



107468, 



ri ■ = 1796678230, ^ = 413605561988912, 



.(2) 



1768256505302499935380 . 



(14) 



We are not able to recognize this fast growing sequence. 
Increasing a, the growth is even faster. 
The numerical values of Rn^ : 

Ri = 0.5, i?2 = 0.507813, R3 = 0.409958, (15) 
i?4 = 0.418322, i?5 = 0.367356, Re = 0.374443. 

clearly display even-odd oscillations that prevent us from 
any precise scaling analysis as the previous one for the 
entropy. Multiplying this by n^^^ (see Eq. ((T^ ) we get 

0.5, 0.604, 0.540, 0.592, 0.549, 0.586, (16) 
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FIG. 2: Spectrum of p„ in the TD limit, the dots represent 
the first and last element Eq. (|6]). Note that these dots sit at 
3/4 (2/3) of the spectrum for odd (even) n. 



that tend to approach a constant value confirming the 
CFT scaling. A systematic study of these oscillations 
requires the analysis of larger values of n, that are not 
accessible to the present method and for which we per- 
formed DMRG calculations that appear elsewhere 14 1. 
For periodic boundary conditions, these oscillations are 
only present for a 7^ 1, and not for the entanglement 
entropy. Thus they have a different origin than those 
found elsewhere for Sn that are due to boundary effects 
The same kind of oscillations have been found 
numerically for the multi-interval entanglement [21]. It 
is rather natural that all the elements of Pn{L) oscillate 
as a consequence of the tendency to antiferromagnetic 
order of the XXZ chain at A = —1/2, and consequently 
most of the averages that are calculated from them are 
expected to oscillate as the moments do. Why and how 
these oscillations cancel between each other only for the 
von Neumann entropy is still mysterious. 

i?!"-* is not the only scaling quantity that can be rep- 
resented as sequence of rationals. A good alternative is 
given by the central values Qr"'' = Rl"\2n + I), which 
grows more slowly. In fact, in the matrix Pni'in + 1) 
the common denominator is the inverse of E-(2n + l,n) 

in Eq. As a typical example we report Qn' — 

qnE'i{2n+ l,n): 



qi =5, q2 = 327, qs = 159502, = 680263760, 
95 = 22821555833635, qe = 6408136183930928388 .(17) 

(2) 

Unfortunately, although it grows slower than r„ , we are 
also unable to recognize the sequence, but the guessing 
seems less prohibitive. 

The spectrum of p„ in the TD limit is shown in Fig. [2l 
Several interesting features are evident in the plot. For 
odd n all the eigenvalues are doubly degenerate, while 
for even n only some. The spectrum at n is roughly 
repeated at the bottom and at the top of the spectrum 
at n-f 2 with some new structure in between. This is also 



0.15 
0.12 
0.09 
0.06 





o ZZ coupling X 2 
□ XY coupling 


8 

I.I. 


B 



0.2 



0.4 0.6 



0.8 



FIG. 3: Dependence of the coupling constants of the entan- 
glement Hamiltonian on the position in the block. Squares 
(circles) refer to the XY (2 times ZZ) couplings, while the 
full lines are only guides for the eyes to connect points at the 
same n. 
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what happens for free fermions jj (A = in Eq. ([3])). 
Thus the interaction appears to change only quantitative 
features of the spectrum and not the qualitative ones. 
This is highly non-trivial because non-zero A introduces 
well-known strong non-perturbative effects. The smallest 
eigenvalue scales like const" (i.e. the top of Fig. [5] is 
almost a parabola). Such scaling is known to be true for 
the "all up" eigenvalue 1], i.e. to the EFP that 
is marked as a dot in the figure. 

The logarithm of the density matrix can be interpreted 
as an effective Hamiltonian for the subsystem through 
the natural definition Pn{L) = e~^"'^^^ at a fictitious 
temperature T = 1 (as observed independently Q)- This 
effective Hamiltonian can be written down exactly only 
for free fermions 0, but it is extremely difficult to ob- 
tain even part of it for interacting systems, because it 
requires the knowledge of the full density matrix. The 
present study gives this unique opportunity (it is enough 
to diagonalize the density matrix, taking the logarithm 
of the eigenvalues and then use the diagonalizing trans- 
formation to go back to the original spin basis). In the 
TD limit we find that the largest terms are a diagonal 
interaction, JjSjSj^^, and a nearest neighbor exchange 
JjSj' S~^_j^+hc with a ratio —0.55 that is almost the 
same as in the original model —0.5. All other terms (ex- 
changes over larger distance, exchange of more than one 
pair of spins, multispin interaction) are more than one 
order of magnitude smaller. The couplings depend on 
the position roughly quadratically 

j;(„)^A^^^, (18) 

with A ^ 0.6. The parabolic dependence of the coupling 
constants is shown in Fig. [3) 

The symmetrized density matrix. The degeneracy of 
the ground-state at finite L gives another unique oppor- 
tunity: the symmetrized density matrix 

p^ = l(|*+)(*+H-|*_)(*_|), (19) 



FIG. 4: Entanglement entropy of the symmetrized density 
matrix S^{L) versus n/L and its asymptotic behavior. The 
full line is (log n) /3 -1- c'l that is an effective description up to 
n/L ~ 0.5. The dashed line is the finite-size CFT prediction 
for pure state (logL/7rsin7rn/Z/)/3 -I- c'l that clearly does not 
work for x > 0.2. The dashed-dotted curve is the heuristic 
formula Sli'^L = n/x) ~ 1/3 log [n(l - ylsinh*a;)] with A ~ 
0.48 that works well for all values of n (except very small n 
and L — n). Inset: Scaling of the mutual information M„ 
as function of n/L compared with the symmetrized heuristic 
guess. 

corresponds to minimum energy and it is a zero temper- 
ature mixed state, has no interpretation in CFT, and 
so it is a new quantity. In the TD limit S{p^) = S{pn). 
In Fig. [4]we report the exact results for the entanglement 
entropy S^{L) of this mixed state. Up to n/L ~ 0.5 they 
are well described by S^{L) — (logn)/3 + c'^, {indepen- 
dent of L and with the same c[ as before) as shown by 
the full line in the plot. Note that a finite size correc- 
tion in the form of a sine (dashed-line) does not work 
for X > 0.2. However, a good collapse is observed for all 
X. (the last points that do not fall on the master curve 
correspond to rt = L — 1 and are not expected to scale) . 
We find that the corrections to the log n behavior are of 
the order of {n/L)'^. The collapsed data are fit remark- 
ably weU by si'\L = n/x) ~ 1/3 log [n(l - ^sinh^'a;)] 
with A ~ 0.48. We do not claim that this form has any 
justification. 

In a mixed state the entropy ([1]) is not a good measure 
of entanglement because it mixes classical and quantum 
correlations. From quantum information we know that 
an appropriate measure is the mutual information 

M„ = Sn + Sl-u - Sl . (20) 

The calculation of Af„ is more difficult because even for 
small n we need Sl-u that can correspond to a large 
block. The data available from the exact ground state 
(up to n = 13, L = 21 plotted in the inset of Fig. 3]) 
collapse and define a universal function that is described 
by the symmetrized version of the heuristic Sn \L) in- 
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troduced before, plotted as a full line in Fig. ID 

Such mixed zero-temperature states are present ev- 
ery time that the ground-state is degenerate at finite 



L. Among these, supersymmetric lattice models [22| 
are very interesting and they could be understood more 
deeply in this framework. 

Conclusions and perspective. We presented explicit an- 
alytic expressions for the reduced density matrix of the 
XXZ spin-chain at A —1/2. From these matrices we 
built several sequences of integers that encode the scal- 
ing behavior. From the exact density matrix, we con- 
structed the entanglement Hamiltonian, a result that is 
not easily accessible to other approaches. We found the 
remarkable property that this Hamiltonian is dominated 
by nearest-neighbor terms of the same form as the orig- 
inal Hamiltonian. Finally, because the ground-state is 
doubly degenerate, we could study the entanglement of 
a zero-temperature mixed state, showing a very different 
behavior from a pure one. 

These results are unique in two ways: they concern 
the full density matrix rather than only its entropy, and 
they are completely exact rather than asymptotic. They 
suggest many questions that can be posed also numeri- 
cally for more general systems. They show explicitly that 
even for modest block size, asymptotic behavior is evi- 
dent, eventually allowing for the study of further finite 
size corrections. Finally, since the ground state is known 
to have strong connections with combinatorial problems 
of current interest d, [l^, [l^ [l^ , it is expected that fur- 
ther properties of the density matrix for arbitrarily large 
blocks can be found by combinatorial methods. 
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